Nflation: observable predictions from the random matrix mass spectrum 
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We carry out numerical investigations of the perturbations in Nflation models where the mass 
spectrum is generated by random matrix theory. The tensor-to-scalar ratio and non-gaussianity 
are already known to take the single-field values, and so the density perturbation spectral index 
is the main parameter of interest. We study several types of random field initial conditions, and 
compute the spectral index as a function of mass spectrum parameters. Comparison with microwave 
anisotropy data from the Wilkinson Microwave Anisotropy Probe shows that the model is currently 
viable in the majority of its parameter space. 



PACS numbers: 98.80.Cq 



I. INTRODUCTION 



The Nflation model of Dimopoulos et al. [H corre- 
sponds to a collection of uncoupled massive fields which 
drive inflation via the assisted inflation mechanism 
The existence of multiple fields is motivated by the ax- 
ions of string theory, and their presence enables sufficient 
inflation to be obtained without using super-Planckian 
field values. Models of this type had first been consid- 
ered by Kanti and Olive ^ and then Kaloper and Liddle 
[H in the context of Kaluza-Klein models. They showed 
the massive fields evolve faster to the minima of their own 
potential and light ones later. Easther and McAllister Q 
introduced random matrix theory as a way of computing 
the possible distribution for masses in the Nflation model. 
For some related constructions see Ref. @. 

It is obviously important to develop observational pre- 
dictions from such models. This has been thought diffi- 
cult, because in multi-field models the predictions depend 
in general upon the field initial conditions as well as the 
model parameters. Easther and McAllister 0] only stud- 
ied two types of initial conditions, where either the field 
values or the field energy densities were equal. Neither is 
well motivated physically. However, in the case of an ex- 
ponential mass spectrum, Kim and Liddle [7] showed that 
provided there were enough fields, with randomly-chosen 
initial conditions, the observational predictions become 
essentially independent of initial conditions again. The 
reason is that with enough fields, the space of possible 
initial conditions is well sampled by a single realization, 
and they named this the 'thermodynamic regime'. 

The purpose of this article is to apply the random ini- 
tial conditions approach of Ref. [?! to the random matrix 
mass spectrum of Ref. , in order to explore its obser- 
vational predictions and test its viability. 

The main observables are the density perturbation 
spectral index no,, the tensor-to-scalar ratio r, and the 
non-gaussianity parameter /nl- Some quite general re- 
sults are already known, applying to arbitrary mass spec- 
tra and initial conditions provided sufficient inflation is 
obtained. Alabidi and Lyth [§| showed that r always has 
the same value as in the single-field case, and Kim and 
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Liddle j^] showed that the same was true of /nl (this 
having previously been shown for two fields in Ref. (1C| ) . 
Only ns has model and initial conditions dependence, 
and it has been shown that its predicted value cannot be 
larger than the single-field value [111, There are some 
further generalisations of these results 1, [l^] . 



II. THE RANDOM MATRIX THEORY 

Once one considers the scalar fields in inflation as ax- 
ions in string theory, their masses can be written in ma- 
trix form, which depends on specific details on compact- 
ification (following Ref. we also assume that higher- 
order terms can be neglected, so that what are really 
cosine functions can be approximated as massive uncou- 
pled fields) . The shape of the mass distribution depends 
only on the basic structure of the mass matrix, which 
is specified by the supergravity potential. In the sim- 
plest assumption, the entries in the mass spectrum are 
independent and identically distributed, i.e. a random 
matrix. The fields can be uncoupled by diagonalization 
of this matrix, with the mass spectrum given by the dis- 
tribution of eigenvalues. The distribution of the eigen- 
values for random matrices of this kind is characterized 
by the Marcenko-Pastur law J^] when the matrices are 
large. The distribution function depends on a parameter 
/?, the ratio of the number of axions to the dimension of 
the moduli space. Easther and McAllister 0] devised this 
formalism and computed the observational predictions in 
terms of /3 for specific choices of initial conditions for the 
fields where the field values or the energy densities are 
identical (see also Ref. (T^V 

We follow their notation for the mass spectrum, and 
label the average value of the mass-squared {m?) = fiT? . 
We will throughout put to = 10~^ Mpi where Mpi is the 
reduced Planck mass; a simple rescaling of to would be 
sufficient to match the observed normalization of per- 
turbations. The mass spectrum is determined by two 
quantities, the number of fields A^f and a parameter (3 
which governs its shape. The shape parameter (3 lies in 
the range zero to one, with values around 1/2 perhaps 
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FIG. 1: The mass-squared distribution for 1000 fields; the 
dashed-dotted (black) flat line is for /3 = 0, the dotted (blue) 
for /3 = 0.3, the solid (red) for /3 = 0.5 and the dashed (green) 
for /? = 0.9. 
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FIG. 2: The spectral index ns (from top to bottom, /9=0, 
0.1, 0.3, 0.5, 0.7, 0.8, 0.9 and 0.95). The dotted line shows 
the observational lower limit on no, from WMAP3. 



the most plausible [1, 5]. Figure [T] shows the shape of the 
spectrum for A^f ~ 1000 and some values of /3. 



III. OBSERVATIONAL PREDICTIONS 

For a set of uncoupled fields with quadratic potentials, 
the number of e-foldings A^, in the slow-roll approxima- 
tion, is pl| 



TV: 



(1) 



where is the z-th field. For our initial results, we 
choose the field initial values randomly from a uniform 
distribution in the range to A/pi. The total number 
of e-foldings is then accurately given by the linear rela- 
tion TVtot — A'f/12 this result holding for an arbitrary 
mass spectrum. We will throughout assume the observ- 
able scales crossed outside the horizon 50 e-foldings be- 
fore the end of inflation. Accordingly, sufficient inflation 
requires a minimum of around 700 fields. 

In passing we note that increasing the number of fields 
increases the energy scale at the end of infiation, and 
hence might increase the number of e-foldings relevant 
to observable perturbations. However this effect is well 
within current uncertainties from the unknown behaviour 
of the Universe between inflation and nucleosynthesis. 

Wc follow the usual formulae for the observational pre- 
dictions of Vn^ "s, and /nl- These are [1, SB 
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where nii and are the i-th mass and the sec- 

ond term of the nonlinearity parameter /nl respectively 
(the first contribution to /nl is model-independent and 
small [l^ [i^). More detailed calculations of the non- 
gaussianity have confirmed that this is the leading con- 
tribution to the non-gaussianity [l3] ■ We have simulated 
these observables and will show in particular rig to study 
the effect from this mass distribution. 



A. The spectral index 

The most important observable in Nflation models 
is the spectral index, which we calculated as described 
above and show in Fig. [2l We see that it depends on 
both model parameters A'f and (3. 

The curves have a generic shape where they first dip, 
and then, after a minimum typically around N{ — 1500, 
increase to join the single-field value. The single-field 
value is always obtained in the case /3 = 0, in which all 
masses become identical and this result is well known 



(e.g. Ref. [Ujl. These curves look a little different from 
Fig. 4 in Ref. [3] , studying an exponential mass spectrum 
where the values asymptoted to constants for large Ni. 
This is just because iVf is defined here in a rather different 
way; in Ref. [7] increasing it added new fields at the top 
of the mass range, and if they were heavy enough they fell 
to their minima before observable scales left the horizon. 
The distinction here is that when we increase N{ we are 
packing more fields into the same mass interval, forcing 
us to the equal-mass case. That the curves begin to rise 
after iVf ~ 1500 indicates that fields heavier than the 
1500th have typically reached their vacuum state before 
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FIG. 3: As Fig. (2] but now showing sets of curves for two 
choices of initial condition distributions. The dashed lines re- 
produce the curves from Fig. [2] and the solid ones show initial 
conditions with (p^ /Mpj drawn from a uniform distribution. 



the 50 e- foldings point. 

To compare with observations, the analysis in Ref. (l6| 
shows that for r = 8/50, the 95% confidence lower limit 
for ng lies at about 0.93, shown as the dotted line in the 
figure. Provided (3 < 0.8, the spectral index is always 
large enough regardless of the number of fields. For large 
/3, the model is excluded only for a range of Nf running 
from about one thousand to several thousand. Accord- 
ingly, most of the model parameter space is currently 
viable. 

The results shown in Fig. [5] are the mean values over 
realizations of the initial conditions. Additionally we find 
that the spread in ng values is quite small; the standard 
deviation of ng is never more than one percent of its 
displacement from unity. This confirms the existence of 
the thermodynamic regime for this mass spectrum, just 
as in Ref. Q. 

Although the above analysis shows that for our chosen 
initial distribution, the observational predictions are in- 
dependent of realization, one might further ask whether 
there is dependence on the choice of that distribution. 
In order to test that, we carried out two further series 
of simulations. In the first series, we took (/)f/M|i to 
be chosen uniformly between zero and one. The results 
are shown in Fig. [3] The curves are shifted to a smaller 
number of fields, retaining both their shape and minima. 
We found that when the number of fields is multiplied 
by 3/2, making the total number of e-foldings the same, 
these solid curves become matched to the dotted ones. 
This can be related to the different initial distributions 
by an approximate analytic argument given in the Ap- 
pendix. 

In the second series, we chose the field energy densi- 
ties to be uniformly distributed, with the lightest field 
ranging from to Mp\ and the other fields given smaller 
ranges to ensure the same expected energy density. This 
is distinct from Ref. Q who gave each field the same 
energy density, rather than the same only on average. 
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FIG. 4: As Fig. O but showing results for uniformly-chosen 
field energy densities. Only f3 = 0.3, 0.5, 0.7, 0.8 and 0.9 
(from top to bottom) are shown. The left-hand end of the 
curves corresponds to models only just achieving 50 e-foldings 
in total. 

In this case the more massive fields start closer to their 
minima, and a greater total number of fields is needed to 
get sufficient inflation. The results are shown in Fig. 31 
The curves have a rather different shape, and there is 
some difficulty in achieving a minimum of 50 e-foldings 
unless N{ is large. The analytic analysis given in the 
Appendix does not apply in this the field ini- 

tial values are correlated to the masses. Nevertheless the 
overall conclusion is the same: that only for large /3 is 
there a danger that the models are ruled out through 
predicting too small a value for ng, and even then only 
in a narrow range of A^f. The expected value (3 ~ 0.5 is 
therefore viable for all initial condition distributions we 
tested. 



B. Other observables 

Concerning the tensor-to-scalar ratio r and the non- 
gaussianity parameter /nl, unsurprisingly the results are 
same as in our previous work 0iH ■ Comparing to the an- 
alytical results of Eqs. (O, we confirm numerically 
that r is indeed independent of the model parameters, 
r = 0.16, and that /nl — 2/A^ which will never be suf- 
ficiently large to be detected. All simulations of r and 
/nl individually give those results, without needing to 
average over realizations of the initial conditions. 



IV. CONCLUSIONS 

We have carried out a detailed numerical investigation 
of the inflationary perturbations in Nflation models with 
the mass spectrum of random matrix theory. We pro- 
duced the predictions for the spectral index as a function 
of both the number of fields N{ and the distribution pa- 
rameter (3 in the mass spectrum. We have found that 
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the model remains viable in the majority of its parame- 
ter space, and that the thermodynamic regime, where the 
predictions become independent of the initial condition 
realization, holds for this spectrum as for the exponential 
case [3. 
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APPENDIX: ANALYTIC EVOLUTION 

This appendix is based on unpublished work by 
Richard Easther, whom we thank for providing details. 

In the main text we considered different initial distri- 
butions for the scalar fields. When the initial values of 
(j)i/Mpi and c/yf/Mp^ are randomly chosen uniformly be- 
tween zero and one, the expectation value a = 
is 1/3 and 1/2 respectively (other choices will give a dif- 
ferent a). Eq. ID) becomes A'tot — N{a/A. This difference 
by a factor 2/3 is closely related to the shift in the curves 
shown in Fig. [3) 

This can be shown via an analytic approximation to 
the evolution, in order to extract estimates of the sum- 



mations 



and X;f*' 



which appear in the 



expression for the spectral index. Using the slow-roll ap- 
proximation, the field value is 



0,(i) = 0,(io)r™'/''(O, 



(7) 



where T{t) and b have the same definition as in Ref. Q: 
they are the ratio of the value of the heaviest field at 
the time t to its initial value, and the upper limit of the 
probability distribution of the mass-squared spectrum, 
respectively. Also is the initial time. Then using the 
exponential function with c — c(t) = 21n[T(t)]/&, we can 
rewrite the summation terms such as 



rn^ifl^iit) = X] ™?'/'?(*o) exp[m^c] 



(8) 



If one ignores correlations between the mass dis- 
tribution and the initial field distribution, so that 
(mf exp[rn^c]0f (to)) = {mf ex.p[mf c]) {(f>'f{to)), then with 
the Marcenko-Pastur distribution, the expectation value 
is given by 



where 2^1 is the hypergeometric function. This follows 
from Eq. (6.14) of Ref. [5|, where the summation can be 
rewritten as a hypergeometric function. 

Therefore the summation terms with the expectation 
values of field initial conditions and of distribution of 
mass spectrum are 
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= 47Vtot Af^^m^ J2 2Fii~i, - 1; 2; /3) I , 
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00 I 

= ANt^^Ml.m^ m^' 2^1 (-z - 1, - 2; 2; /3)^ . 

1=0 

Note that t is related to Ntot—N; that is, t counts forward 
from when the fields start to evolve from their own initial 
conditions, while N is the number of e- foldings before the 
end of infiation. 

The ng data in Fig. [3] are taken at different times, but 
at the same N. Using Eqs. PH)) and pT|) . the ng formula 
becomes 



, 1 f(t,P) 



where 
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Here the function f{t,l3) represents how rig evolves with 
respect to time (or TV). Even though it is hard to evaluate 
it in general, it will have the same value whenever t and 
/3 are same. Hence in cases with same A'tot, but not 
necessarily the same N{ or a, then ng must be same. 
This explains why the solid curves in Fig. [3] are matched 
with dotted ones after multiplying 3/2 to Nf, because 
this makes A'tot the same. Different a shifts the ng curves 
along the A'f axis, with bigger a shifting the curves in the 
smaller A^f direction. 

We can conclude then that, insofar as the approxima- 
tions hold, (3 < 0.8 should satisfy the current observa- 
tions regardless of the field initial condition distribution, 
while for larger (3 some choices of N[ will be ruled out 
and that those Nf values depend on the a value of the 
initial condition distribution. However the above analy- 
sis relies on using a slow-roll approximation for all fields, 
which is likely to become increasingly inaccurate as in- 
creasing numbers of fields evolve towards their minima, 
and assumes the mass and initial field value distributions 
to be uncorrelated. The latter approximation fails badly 
for the uniform field energy density initial conditions. 
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